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Abstract 

We consider the Cauchy problem for the weakly dissipative wave 
equation 

Dv+-^— v t = 0, t>0, 
1 ~\~ t 

parameterized by \i > 0, and prove a representation theorem for its 
solutions using the theory of special functions. 

This representation is used to obtain L p -L q estimates for the solu- 
tion and for the energy operator corresponding to this Cauchy problem. 

Especially for the energy estimate we determine the part of the 
phase space which is responsible for the decay rate. It will be shown 
that the situation depends strongly on the value of fi and that /i = 2 
is critical. 
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1 Introduction 

We are interested in properties of solutions to the Cauchy problem 

av + -^- t v t = 0, v(0,-)=v 1 , v t (0,-)=v 2 (1.1) 

to data vi,v 2 6 S(R n ). 
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I. On the one hand this problem is the model case of weak dissipation. 
The problem is investigated by this reason, but not well understood. This 
shall be underlined by one result. If we consider the usual hyperbolic energy 
of a solution 

E(v,t) = \ [ {v 2 + (Vv) 2 )dx, (1.2) 

z JR™ 

latTfj ], [Ues7S| ] that 

E(v;t) = 0(t~ a ), a = min{2,/x}. (1.3) 

The method of proof are weighted energy inequalities. They provide a pow- 
erful tool of proving such a-priori estimates without knowing details on the 
solution representation. 

So far it is not motivated why this change in the decay order for \x = 2 
appears. We will give an alternative proof of such a-priori estimates using an 
explicit representation of the solution operator by means of special functions. 
Later on we will see that different techniques are necessary to understand 
small and large values of the parameter fi. 



it is known from 



II. On the other hand this equation is related to wave equations with 
increasing in time propagation speed. If we are interested in solutions to 
the Cauchy problem for the operator 



A? - A 2 (r)A, 



(1.4) 



investigated e.g. in | RY98H , we can apply a change of coordinates to reduce 
it to an equation with constant speed of propagation and dissipative term. 
We introduce the new time variable 



such that d T = X(r)dt and d. 
the operator rewrites as 



t 

TO 

2 _ \2/V\a2 



X 2 (r) 



X(s)ds, 

\ 2 {t)8 2 + \'{T)d t . Thus in the variable t 
\'(r) 



□ + 



A 2 (r) 



dt 



(1.5) 



Especially for the coefficient A(i) = (1+i) , investigated by M. Reissig in 
[ Rei97ll , we get with r = (1+f) W -1 the relation t = 1/(^+1) (l+r) m -l 
and therefore the equivalent problem 

£ 1 



□ + 



£+ll+t 



dt, 
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i.e. Q with fi = £/(£ + 1) G (0, 1). 

Similarly for X(t) = e t we get the relation to our equation with parameter 
fi = 1. This Cauchy problem was recently studied by A. Galstian in (GalOOj. 



This paper extents the results of [ Rei97| and [|Gal0C|1 and provides in the 



case of the latter one even sharper estimates. 



In all the papers [Rei97], [GalOO] and | TT80 | the Cauchy problem was 
reduced to the confluent hyper geometric equation by change of coordinates 
in frequency space. We will follow this idea. 

III. It should be remarked that the related singular problem in one space 
dimension 

vtt + ~vt - v xx = 

is related to the free wave equation for radial data (understanding x as time 
variable and t as radius). This follows directly from the representation of 
the Laplacion A in polar coordinates 



,.n n — 1 _ 1 , 
A = d 2 r + d r + —A 



s, 



where As is the Laplace-Beltrami operator on the unit sphere S 71 " 1 of R n . 

Thus this special problem generalizes the Euler-Darboux equation to 
non- integral /i. 

In |Wei54] a solution representation for the multi-dimensional analogon 



k 

vtt + -v t - Av = 

was obtained using recursion formulae and a generalized method of descent. 
The difference to our representation is that we use Fourier multiplier and 
work in frequency space while A. Weinstein used integral formulae in phys- 
ical space. 

IV. In the scope of this paper are two main points. 

At first we give an explicit representation of the solutions to the Cauchy 
problem (LI) in terms of Bessel functions. This will be done in Section ||| 



Afterwards we use this representation to describe the Fourier multiplier 
corresponding to 

(yi,v 2 ) i-> (Vv,v t ) 

in order to derive energy estimates in the (strong) sense of a norm estimate 
for this operator. The same is done for the solution operator of the Cauchy 
problem itself. This is the main contents of Section ||. 
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In Section || we give a further a-priori estimate including a non-trivial 
source term. 



2 Multiplier Representation 

Reduction to Bessel's equation. At first we construct the fundamental 
solution of the corresponding ordinary differential equation in the Fourier 
image. Let v(t, £) = ^-^[v], 

t)(t,0 = (27r)-5 / e- ix <v(t,x)dx. (2.1) 



Then {) satisfies the ordinary differential equation 

vtt + \^\ 2 v + jr^vt = 0. (2.2) 



Following [TT80], we use the relation of this differential equation to Bessel's 
equation to construct a system of linearly independent solutions. 
We substitute r = (1 + and obtain 

d dr d . d . 

= did? = I«ld7 < 2 ' 3 »> 

d 2 d 2 

and hence the differential equation 

|£| 2 ^ + 1^ + 1^1^ = 0. 

After multiplication with |£|~ 2 we get the |£ [-independent differential equa- 
tion 

v TT + £y T + y = 0. (2.4) 

If we make the ansatz v = t p w(t) this leads to 
=p(p - 1)t p ~ 2 w + 2 P T p - l w' + t p w" 

+ - {pT P ~ l W + T P Vj') + T P W 
T 

+t p ~ 2 (r V + (/i + 2p)ru/ + (r 2 + p(p - 1 + fi))w) , 
i.e. by the choice of p + 2p = 1, 

P=-^, (2.5) 
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and hence p — 1 + \i = —p we get Bessel's differential equation 

r V + tw' + (r 2 - p 2 )w = (2.6) 



of order ±p. A system of linearly independent solutions of fl2.6| ) is given by 
the pair of Hankel functions 7V^(t). For details we refer to the treatment in 
[WaE§]. 



Hence 

w + {t) = tPH + p (t), w-(r) = t"H-{t\ (2.7) 

with p determined by ( p.5| ) gives a pair of linearly independent solutions of 
(p.4|) and after performing the substitution r = (1 + of (|2.2|). 



Representation of the Fourier multiplier. We are interested in a pair 
of special solutions *i(£, to;£), ^2(*)*o>£) °f Q2.2|) to initial conditions 



* 1 (t 0) <o,0 = l, $*i(*o,*b,£) = 0, (2.8a) 

$2 (to, to, = 0, d t $2(t , t , = 1, (2.8b) 

where the parameter to > — 1 describes the time level, where the initial data 
are given. We collect these $j in the matrix 

$(Mo ' e) -U*i(Mo,e) ^ 2 (Mo,oJ- (2 ' 9) 

Then for the vector (v(t, -),vt(t, -)) T we obtain the relation 

$3) =•"■*••> (M- (2 ' 10) 



The matrix $ is called the fundamental matrix of ( |2.2[) . 
For iu±(i, £) we have the following initial values 

w + (t ,0 = (1 + t ) p |£| p ^((l + t )|£|), (2.11a) 

^+(t , o = (i + toner 1 ^-! ((i + to) ici) , (2.nb) 

w.(t ,o = (i + t ) p iei p ^;((i + to)iei), (2.iic) 

d t w^t ,0 = a + to) p \Z\ p+1 n;_ .{(i + toM), (2.iid) 
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which follow straightforward from the recursion formulae for Bessel functions 
| Wat22| , §3.6]. For instance we have 

^[(i+w^((i+*)iei)] U 
= po. + ty-^H+tti + m\)\ t=to + (i + ty\tv\H + p )>{{i + m\m\ t=t0 

= a + foner 1 (jf|tt?((i + *o)iei) + cw+)'((i + mi) 

= (i + t r\^ +1 n;_ 1 ((i + toM\). 

From these initial values we determine the constants Ci±(to,£) in 
^ i (t,t ,0 = C i+ (t ,Ow + (t,0 + a.(t ^)w.(t,C), z = 1,2, (2.12) 
such that (|2.8|) holds. This means we have to satisfy 



( w+(to,t) w-(t ,£) \ (C 1+ (t Q ,0 C 2+ (t ,0\ = T r91 oN 
\dtw + (to,0 dtw-(to,0) VCi_(t ,0 c 2 _(t ,0y ' 



Hence we have 



/Ci+(t ,0 C 2 +(io,£A = 1 / ^-(to,C) -«>-(*>, 

VCi_(t ,0 C 2 -(t ,0) detW{t ) \-d t w + (t ,0 w+(to,0 , 

(2.14) 



where 
det W(i ) = det 



/w+((i+t )iei) W((i+<o)iei)V 
" ( + oj ICI K-((i+i )iei) (H,-)'((i+t )iei)J 

4? 

= leP'Cl + to) 2 "- 1 , (2.15) 

7T 

because the Wronskian VV of the Hankel functions satisfies 

W(H P + W,H P -W) = de t (^ (2.16) 



sec 



Wat22| , §3.63]. 



Hence we obtain for the fundamental solution 



$i(t, t , o =ci+(*o, ew*, e) + ci-(*o, o«-(t. o 



ficr^a+to) 1 - 2 " 

(1 + io)lCr + X"-i ((! + ( X + OIA"*? ((! + OKI) 

- (i + t y\w +i n+ ((l + <o)iei) (i + tywn- ((i + m\) 



-H+((l + i )KI)^((l + i)|e|) 



(2.17a) 



and similar 

$2(M0,£) 



|7T^i{^-((i+*o)iei)^((i + t)iei) 
-H+((i+< )iei)wj((i+t)iei)). 



(2.17b) 



We collect the results in the following theorem. 



Theorem 2.1. Assume v = v(t,x) solves Cauchy problem fll.ip . T/ien i/ie 
Fourier transform v(t, £) satisfies 

.7=1,2 



where the multipliers are oiwen 6y 



$i(Mo,6 



Z7T , j_, {l+t)P 



and 



$ 2 (Mo,£) 



4" I5I (1 +i )P-! 
in (l+t)P 

~T(i + t )^- 1 



w;-i((i+*o)iei) Wp((i+t)iei) 
wti((i+«o)iei) w?((i+t)iei) 
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Time derivatives. Next we need the time derivatives of these functions. 
Derivation with respect to t leads to 



id 



w+i((i+*o)iei) n+{{i+t)\i\) 
w^i((i+<o)iei) (Wp)'((i+*)iei) 

^((l+^iei) ^((i+^ici) 
wf.i((i+<d)iei) wf.i((i+*)iei) 

using pTC^(z) + z{TL~^)'{z) = zHp_ 1 (z) and similarly for the second multi 
plier. 

Corollary 2.2. T/ie iime derivatives of the multipliers <3?j are given by 



^ 4 l?l (l + i )p-i 



9t*i(t,to,0 
and 



m 2 (l + t)P 
4 l?l (1 + io)"" 1 



4 l?l (l + i )p-i 



w;_i((i+*o)iei) w;_i(a+*)iei) 

w3Li(a+*d)iei) w£-i((i+*)iei) 

w-((i+to)iei) w^i((i+t)iei) 

w P + ((i+to)iei) w+i((i+t)iei) 



It is possible to obtain a similar expression for higher order time deriva- 
tives in the same way. 



Representation by real-valued functions. If we use the definition of 
TL^ by the real-valued Bessel and Weber functions [Wat 2 2 , §3.6] 

nf(z) = j p (z)±iy p (z), (2.i8) 

we obtain an alternative characterization of <1> by real-valued functions 

j p -i«i+t )\(\) jMi+tm) 
^((l+^iei) y P (o-+t)\t\) 



$i(Mo,0 



2 ICI (l + t )P-i 



*a(Mo,0 



vr (i + ty 



2(1 + to)"" 1 



(2.19a) 
(2.19b) 



y,2 a+*) p 

2 ig (1 + io)"" 1 



2 l?l (l + t )P-i 



^>((i+«o)iii) M(i+m\) 
y P ((i+t )\z\) y P ((i+t)\z\) 

^((l+^iei) j P -i((i+m\) 
y P -i((i+to)iei) yp-i((i+t)iei) 

(2.19c) 

Jp((l + to)lel) Jp-l((l + t)|^|) 

^((i+to)lci) 



(2.19d) 
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The representation simplifies in the case of non-integral p using 

y p {z) = cot (pn)J p (z) - esc (pn)J- p (z) 
and the periodicity properties of the trigonometric functions to 

(i + ty 



$i(Mo,0 



$ 2 (Mo,0 



&*i(Mo,0 



$$2(M0,O 



7T 



2 c S c(pvr)|C| II + ^ )p _ 1 
JLp+i^l + to)^) JL p ((l+i)|$|) 

-jp-i((i+t )iei) jp((i+i)iei) 

l csc( ^ (i ( +t )^ 
j_ p ((i+t )i?i) j_p((i+*)iei) 

^((l + *o)l£l) «7 P ((l + t)|£|) 
•I CSc( ^ )le|2 (l+! )p-i 

j_ p+1 ((i+i )iei) j_p+i((i+t)iei) 
^((l+fo^i) jp-^i+m) 



(2.20) 



(2.21a) 



(2.21b) 



(2.21c) 



(2.21d) 



In the first and in the last formula we used esc (pn — n) 



-esc (pn). 



3 Estimates 

We use the isomorphism 

(B) s : L p , s (R n ) L p (M"), 

(D) the pseudo-differential operator with symbol (£) = a/1 + |£| 2 , to char- 
acterize the Sobolev spaces of fractional orderQover L p (R n ), p £ (1, oo). 

The representation of the fundamental matrix <!>(£, 0,£) (as well as the 
knowledge about strictly hyperbolic problems) imply a natural regularity 
difference for the data of one Sobolev order. Therefore we define the follow- 
ing two operators corresponding to the Cauchy problem (|1.1| ). On the one 
hand we are interested in the solution itself. Let 

S(t):(v 1 ,(D)- 1 v 2 )^v(t r ) (3.1) 



^ee e.g. ]AS6l| ] or ]AMS63| for details about Bessel potential spaces. 
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be the solution operator. (Obviously) it maps §(t) : L 2 (M n ,R 2 ) — > L 2 (IR n ). 
On the other hand we are interested in energy estimates. We define the 
energy operator 

E(t) : ((D) Vl ,v 2 ) ^ (d t v(t, •), |D|u(i, •)), (3.2) 

with E(t) : L 2 (IR n ,M 2 ) -> L 2 (M n ,IR 2 ). 

We will give norm estimates for both operators from L p scale to L q with 
dual indices p and q, i.e. p + q = pq. 

Remark 3.1. The norm of the data "contains" (D) v\ and not only the ho- 
mogeneous component |D|«i. The dissipative term has influences like a 
(possibly negative) mass term and thus transfers information of the size of 



v to the energy v', cf. Appendix A. 2 for the transformation of the problem 
to a Klein-Gordon equation. 

Properties of Bessel functions. To obtain norm estimates for the op- 
erator families §(t) and E(t) we have to review some of the main properties 
of Bessel functions for small and large arguments. For details we refer to 
[ Wat22| , §3.52,§10.6 and §7.2]. 



Proposition 3.1. 1. The function 

K(t) = t- v J v {t) 
is entire in v and r. Furthermore, A u (0) ^ 0. 

2. Weber's function y n { T ) satisfies for integral n 

y n (r) = -Mr) log r + A„(r), 

7T 

where T n A n (r) is entire and non-zero for r = 0. 

3. The Hankel functions ri„{r) with r > K can be written as 

7it(r) = e^at(r), 

where a^j € S~z(K, oo) is a classical symbol of order —1/2. 
4- For small arguments, < r < c < \, we have 



\nt(r) 



< 



— log r , v = 0. 
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3.1 L 2 and L p —L q Estimates for a Model Operator 

The model operator. Due to the special structure of the fundamental 
matrix $(f,ioj£) and hence of the multipliers of §(t) and E(i), we consider 
the time dependent model multiplier 



5+1-fc 



(3.3) 



parameterized by A;, s, p,5 £ R. 

Again we can write ^fc a n <j(t, i n terms of the real-valued Bessel func- 
tions of first and second kind. Similar to ( |2.19| ) and ( [2.21 ) we have 



* k ,s, P A t ,0 = m k (O s+1 - k 

= 2i csc(h \a\ k (o s+1 - k 



M\z\) j P+s ui+m\) 
yM\) y P+s ((i+m\) 



(3.4a) 



j- p (\t\) ji p _i((i+t)iei) 
(-l^dei) ^((i+t)iei) ' 

(3.4b) 



the last line holds for p £ Z and p + <5 ^ Z. 

L2-estimates. The structure of the multiplier hints to a decomposition of 
the phase space R+ x R^ into three regions 



Zi 

z 3 



m>K} 

{\£\<K<(l + t)\Z\} 
{(l + t)\t\<K] 



(3.5a) 
(3.5b) 
(3.5c) 



with a fixed parameter K > 0. 
t 



We prove 



\ 22 






Zx 






K \i\ 


^k,0, P ,s(t, 


•) e Loc(M n ) 



(3.6) 
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under certain conditions on k, p and 5 and derive an estimate for 

sup I (0~ S ^k,s, P ,s(t,^)\. 



(3.7) 



Prom the isometry M| ~ (see e.g. | Hor6C , Theorem 1.5], in general 
Mp denotes the space of all Fourier multipliers mapping L p into L q ) this 
gives a norm estimate for the operator corresponding to the multiplier ( |3.3| ) 
as operator H s — > L2. 

Lemma 3.2. It holds ^k,s,p,s{t,-) G L OQ (M n ) for all t if and only if s < 
and k > \8\. Furthermore the estimate 



k,s,p 



< 



f (l+<)-3 )P ^0,|p|-&<-i, 

(1 + ,p^0,H-A;>-i, 
k (l + t)-*log(e + t) ,p = 0,fc<±. 



Proof. We subdivide the proof into three parts corresponding to the three 
zones Z±, Z2 and Z3. 



Z\ We use Proposition 3.1.3 together with the definition of the zone 



K <~]£| < (1 + t)\£\. Thus the multiplier is bounded in Z\ iff s < 0. It 
satisfies 

\^k,s, P ,S^0\ £ (1+*)"* 

un der th is assumption. 



For p ^ we can use |3.1| .4 to conclude 

l*M,p,«(* > oi^[ei fcHp| (i+<r'[er' 



< 



(l + t)"5 ,\p\-k<-l, 
(l + t)M~ k ,| p |-fc>-i. 



For /9 = we have to modify this estimate by the log term 
|*MA^0|^lC| fc log^(l + t)-'|^|-5 



< 



(l+t)"2 ,fc> 

(1 + t)~ k log(e + t) ,k<\. 



2- 
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Z% In this zone we use the representation of ^k,s,p,s{pi i n terms of 



real- valued functions given by (|3,4| ) . For non-integral p and p + 5 we can 
use the representation by Bessel functions of first kind to conclude^] 



k,s,p,S 



(t,0\ 



\t\>j- p (\t\) ((i + m\y p ~ s j P+s ((i + m\) (i + ty +s \z\ k+s 



+ 



p-6\£\k-6 



<(! + *) 



\p\—k 



The condition k > |5| is necessary and sufficient for the boundedness in £. 
For integral values of p or /?+<5 we use Weber's functions and Proposition 



3.1.2. We sketch the estimate if both p and p + 5 are integral. Then we have 



4i, 



k,s,p,8 



io g (i + t)j P mj p+5 ((i + oif i) 



7T 



+ 2t|*|*<fl 



Jp(l^l) Jp +5 ((l+i)l£l) 

A P {\H\) A p+d {(l + m) 



and hence 



I^Wfo 1 S lo s( e + *)(! + *)- |phfc + (1 + i) lphfc 



< 



'(l+t) |p| " 



[(l + t)- fc log(e + t) ,p = 0. 
If only one of both indices is integral we have to mix the representations.^] □ 

L p — L g -estimates. We consider estimates for the model operator 

u(x) h-> [* k ,., p ,s(t, £)«(£)] C«0 (3-8) 
from L P)r (R n ) to L g (]R n ), (p, g) a dual pair and r a sufficiently high regularity. 



By estimating the difference structure of the multiplier by triangle inequality we do 
not loose information. For non-integral p the leading terms of the series expansions do 
not cancel. 

3 For later reference we use only integral values of S. 
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Lemma 3.3. Assume p G (1,2], q such that pq = p + q. Let further k > \8\. 
Then the model operator ( |3.8[) satisfies the norm estimate 



|*fc, a ,P,«(*. D )llw 



< 



n-1 I 1_1 \ _i 

(1 + i) 2 v p »v 2 

( i + t)-»(H)+w- fc 

1(1 + *) 



i-- )+Oe-k,, , , , u1 

p «/ (log^ + t)) 1 



p = 0,<i<±+e,e>0 



/or d 



n+l / 1 
2 U 



+ A; — |p| and r 



77, 



constant in the last case is given by 



1 — 1 

v <? 

n+l 



+ s. T/ie interpolating 



2e+l-2k \p 



Proof. Again we decompose the phase space into three zones. For this we 
use a smooth cut-off function tp € C°°(M+) with ijj' < 0, ip{r) = 1 for r < 1/2 
and ip(r) = for r > 2. Using this function we define 

0i(t,£) = l-V(lfl/#), 

^(t,e) = v(iei/-Kr)(i - ^((i + 
3 (t,e) = V'(i£i/^((i + t)i£i/^), 

such that £) + 02 (i, £) + 03(*> £) = 1 an d supp0i ~ Zj. Thus we can 
decompose the multiplier ^k,s,p,8 into the sum X^i=i 2 3 0*(^> O^fc.s.p,^; £) 
and estimate each of the summands. We prove the estimate for r = 0, i.e. 
we restrict the proof to the corresponding value 



1 1 
-n l 

n q 



From Proposition |3.1| we conclude the representation of 
as sum of two multipliers of the form 



with symbols a E S s+ i(K, 00) and 6 G S~2(K,oo). We follow [|Bre75(| to 
estimate the corresponding Fourier integral operator. The key tool is a 
dyadic decomposition together with Littman's Lemma, [ Bre75| , Lemma 4]. 
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Let X e C, 



o°°( 



be nonnegative with support contained in [1/2, 2] and 

oo 

^ x(2 j r) = 1, for r / 0. (3.9) 



j=-oo 



Such a function exists from |H6r6C, Lemma 2.3]. Let further Xj(£) = 

We obtain an L p -L q estimate for this operator by interpolating L\-L c 
and L2-L2 estimates using [ Bre75 , Lemma 3]. For this we define 



X J (e)e ± ^ l a(|e|)6((l + t)|e|) 



and 



X J (e)e ± ^a(|e|)6((l + t)|e|) 



(3.10a) 
(3.10b) 



and estimate these norms. For all j < we have Ij = Ij = 0. 
For Ij we perform the substitution £ = 2 3 Kr] and obtain 



Ij < C2 ]n 



< C2 jn (l + 2 j Kt)~^ ^ \\D a a(2 j K V )b((l + t)2 j Ki])x{\i 

\a\<M 

< C2 jn (l + 2 j KtY " 



n-1 
2 



sup (2 i K|r ? |) s+ ^-l Q l2 J 'l a l((l + t)2 J K|? ? |)-5-l/ 3 l((l + t)2 J ') 



|a+/3|<A/ 



l/2<|rj|<2 



< C2 J '( n+s )(l + 2 J ia)~ V(l + t)-§ 

by Littmann's Lemma. From C# (1 + *)<(! + 2*1*"*) < C' K 2 j (l + 1) we get 
finally 

Ij < C2^ n+S \l + t)~^ . (3.11) 

For Ij we obtain 



/j<C sup M 2jK r))\^ j Krj\ s+ 2\(l + t)2 j Kr]\ 

»?Gsuppx 



(3.12) 



The estimates ( |3.1lD and ( 3.12 ) correspond to Li-L^ and L2-L2 esti- 
mates for the dyadic component of the model operator (|3.8[). Interpolation 
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leads to 

[Mt,t)xj(®*k,.,p,s(t>Z)*(S)]\\ 9 

< C2 j K^) +s )(l +ty^(v~*)~^\\u\\ p (3.13) 

for all p 6 (1, 2], p + q = pq. Finally, we use |Bre75 , Lemma 2] to conclude 
for n ^| — + s < the estimate 

/ 1 i \ i 



t- 1 [Mt,0^k,s, P At^MO]\\ q < c(i + t) - u *\\ u \\ p . (3.14) 

In this zone we subdivide each summand of the multiplier in two 



Z- 



factors 

\C\ lpl+£ nJ(m 21 (0, (3.15a) 
((1 + i) 1^1) fc_lpl_£ ^ +5 ((1 + ^1) 022 (t , , (3.15b) 

and the remaining constant (l + t)\p\- k+e , where 

foi(ti)= m\/K) and <h2(t,t) = l-iK(l + t)\£\/K) (3.16) 

such that 2 (t,£) = 021 (£) 022 (*,£)• 

The first multiplier is time independent and satisfies 

\^ +£ nj(m 21 (0 ~ ICI |p|+£ ^|(|^|)02i(O 

=(i ± icot iPk)iei 2|p|+ i^r lp| ^ P |(iei)02i(6 

Ticsc|p|7r|^| e |e| l ' ,| J- H (|e|)02i(O, 

=iei 2|p|+£ i^r |p| ^p|(iei)02i(o 
±;iei 2|p|+e log i^i i^r lp| ^ P |(iei)02i(o 
±^ien^i |p| ^ipi(iei)02i(o, ^z, 

where r :fp l l 7±p(r) and t^ p ^A\ p \(t) are entire. By ~ we denote equality up to 
a multiplicative constant here. 

From | Ste7C , §3 Theorem 3] it follows 

ie| £ 02l(|ei) S MP Ve > (3.17a) 
I^| £ log|el02i(!el) €M* Ve>0 (3.17b) 
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for all p € (1, oo). Thus we conclude with the algebra property of multiplier 
spaces that the first multiplier belongs to Mp for p G (1, oo) if e > and 
p^O (or forp = 0ife>0). 

For the second multiplier we prove an L p -L q estimate. For this we use 
again a dyadic decompositon. Let x De like in the discussion of Z\ and 



j>0 



We estimate 



and 



k-\p\-en,± 



xibW 1 + m) k ~ M ~ £ nt s ((i + m\)<i> 22 (t,o 



(3.18a) 
(3.18b) 



(3.19) 



(3.20) 



For j > we have 4>22(t,£,) = 1 on suppxj- Hence using the substitution 
(1 + t)£ = 2 J Krj we get the estimate 



< C2 jn {l + t)~ n 



a(2>Kr,)tl>(\Ti\) 



< C2 jn (l + t)~ n (l + 2^)-V(2^)-i+ fc -|p|- £ 
where a S 5'~2+ fc ~l^l~ e . For In we obtain a similar estimate in the same 



way. 



For Ij we have 



Ij < C sup {2? ' Kr]) k -\ p \- £ -^ 

r;£supp i/j 



(3.21) 



Interpolation leads to 



ut,^+m\) k ' lPh£ nf +s ((i+t)\0H2(t,om 



(3.22) 
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which gives for 

n + 1 (\ 1\ 1 , 



2 



the Lp-L q estimate 



{(i+m\) k - M - E nt s ((i+m\)h2(t,om 



< C{l + t) "Vp *J\\u\\ p . (3.23) 

The 'regularity' e is determined from both multipliers, hence the optimal 
choice is 

e = max {o, ^ (± - ±) - \ + k - \p\} . (3.24) 

under the assumption /> / 0. For p = the choice e = is not possi- 
ble. Therefore for k < ^ we have to exclude that case. We postpone this 
exceptional case. 

Multiplication of the multipliers corresponds to a concatenation of the 
corresponding operators. Hence we have 



q 

i i ' 



F~ [02(t,O*M,P,«5(*.O«(O] 



q 



<G(l + t) M{ -^(F-^-*'- B (F-|)+l''l-*>|| tl || p (3.25) 



for p or p = and A; > 1/2. 



Z2 for p = and k < 1/2 In this exceptional case we get an estimate 
for all dual p and q by interpolation. From 5 + fe > we can follow the 
previously used reasoning to conclude an estimate for dual p, q with 

n+1 (I 1\ 1 



+ /c> 0. 

P ?/ 2 

If we interpolate the corresponding result with the previously proven L2 
estimate in that zone we obtain 

l^- 1 [Mt,S)*k,.Mt,t)W)]\\ q < c(i+ty<^) +9t - k (iog(e+t)) l - s \\u\\ p 

(3.26) 



for 2±i 





-1) 


1 


V 




2 



n + 1 /l 1 



2e + 1 - 2k \p q 
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We use the estimate 



'(l+t)W-* 



;i +t)- fc log(e + t) ,p = 



together with the definition of the zone Z% to conclude the estimate 

< 1 1 03 (*, •)ll 1/ Cl_iyi«llffll*fc ) ».ft*( t ' ^l 00 



< 




(l + t)- fc log(e + i) ,p = 0. 



Under the assumptions on p and g this estimate is weaker than the estimates 
in the zones Z\ and Zi- Q 

Remark 3.2. The zone Z\ determines the necessary regularity of the data, 
large values of |£| are contained in this zone only. This is a difference to the 
reasoning in |Rei97f| , |RY98 | and [ReiOO]. There the pseudo-differential zone 
contains {(t, £)\t = 0}, and hence this zone determines the regularity of the 
data in this case. This is similar to our consideration if we decompose Z\ in 
small t and large t. For large t we can replace ( p.ll ) by 

Ij < C2^ n i 1+S \l + t)~%, 

which gives as conclusion an estimate with less regular data. Hence the 
necessary regularity is determined by the boarder of Z\. 

Remark 3.3. It is possible to extend the estimates to the case p = 1 and 
q = oo by Sobolev embedding. Therefore one has to use slightly more 
regularity, r = n + e. 



Remark 3.4. In Lemma 3J3 we can weaken the assumption k > \6\ in depen- 
dence on the given dual pair (p, q) using the Hardy-Littlewood theorem of 
fractional integration to estimate in Z3. This will not be necessary for our 
considerations. 
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3.2 L 2 and Energy Estimates 

If we compare fl2.17|) with (|3.3|) we obtain the representation 



*i(t,o,o 



/■7T 



(i + t)^i Ap -i,i(t,C), 



Z7T 



— (i +t)"*o,- llP> o(t,e), 



i7T 



^*i(t,o,o = -(i + t)^ 2i i, p _i, (t,C), 



?.7T 



(l+^l.O.p.-i^O 



(3.27a) 
(3.27b) 
(3.27c) 
(3.27d) 



of the entries of the fundamental matrix by our model multiplier. Thus 
we can apply the estimates of Lemma to get a-priori estimates for the 
solution v = v(t,x) of ([Q|). This gives 



K*> Oik ^ 1Mb + \\v2\\h-i < 



((1 + i) 2 " ,pe(0,i), 
log(e + t) ,p = 0, 
1 ,p<0. 



(3.28a) 



For the first derivatives we obtain 

\\v t (t, oi| 2 < ci(i + O^IM!^ + c 2 (i + t ) m ^p-h-n\ 

||Vu(t,-)|| 2 <Ci(l +t) p -3||ui|| H i +C 2 (l + i) max{p "^" 1} | 



V2W2, 



V2W2, 



(3.28b) 
(3.28c) 



which reestablish already known results on the energy decay, lMat77| , |Ues79H 
and |[HN0H , Example 2.1]. 

We collect the above estimates in the following theorem. 

Theorem 3.4. 1. The solution operator §(i) satisfies the norm estimate 



f(l+*) 



TON 



2^2 



< I 



,A*G(0,1), 



log(e + i) ,fi = l, 
1 ,n>l. 



2. The energy operator E(i) satisfies 



\\E(t)\\ 



2 >2 



(l + t)-2 , M e(0,2], 
(1 + i)" 1 ,M>2. 
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Especially for /x = 2 this energy estimate is sharper then known esti- 
mates. For this special case see also Appendix |A.2j . 

Remark 3.5. It should be remarked that it is essential to use the i? 1 -Norm on 
the right hand side of the energy estimate (or the normalization by (D) 1 in 
the definition of E(t), formula fl3.2| )). Otherwise we get for the usual energy 
from Lemma JL2 only the trivial (and in view of this Lemma also sharp!) 



estimate 

E{v-t) <E(v,0). 

This implies that information about the size of the datum v\ is necessary 
for precise a-priori estimates of the energy. 

Remark 3.6. If one is only interested in the decay order, but not the depen- 
dence on the data, one can obtain the slightly better result 

lim(l+t) min ^E(v;t) =0. 



t— >oo 



This result follows from the considerations in [HN01 ] 



Remark 3.7. We see that the decay rates are determined 

• for small /i by the exterior zones Z\ and Z<i and 

• for small fi by the interior zones Z2 and Z3. 

This explains the change in the energy estimate with critical value fi = 2. 

3.3 Lp—L q Estimates 



We apply Lemma 3.3 to the operators S(i) and E(i). Therefore we use the 
representation of the fundamental matrix $ by the model multiplier ^fk,s,p,S 
given in the previous section. 

Theorem 3.5. 1. The solution operator E>(t) satisfies the norm estimate 

\\Ht)\\L p , r ^L q 



< 



,/i6(0,l) 



n-l I 1 11 1 



[1+t) * \p 9) * ,n = l,6> \ 



[I + t y n ^~V +e \\og{e + t)) 1 - 9 >At = l,<5<| + e ,e>0 

(1 +t) 2 \ p q > 2 \p 1J 1 ,^ > 1, 

forp € (1,2], q with pq = p + q, 8 = ^ f± - |) and r = n (J - |V 
The interpolating constant 6 in the third case is given by 9 = ^frj- 
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2. The energy operator E(i) satisfies 

PWIU^l, < (i + t) » u »• u w 

/or p G (1,2], g with pq = p + q and r = n f | — |^ . 

Remark 3.8. The estimates coincide for /i = (i.e. p = 1/2) with the 
well-known L p -L q estimates for the wave equation lStr70fl . For fj, = 2 (i.e. 
p = —1/2) we can reduce the Cauchy problem ( |1 . 1|) to the Cauchy problem 
for the wave equation, see in Appendix |A.2j . 

Remark 3.9. Finally, these estimates are related to L p -L q estimates for 
equations with increasing in time coefficients. By a change of coordinates 
(already mentioned in the introduction) we can transform problem ( |1.1| ) to 

v tt - \ 2 {t)Av = 

with X (t) = (1 + t ) e , i > for (i € (0,1) and X(t) = e* for /i = 1. In 
[ Rei97| ] and [ |GalO0| ] L p -L q estimates for these problems with increasing in 
time coefficients were studied. The estimates obtained here for fi £ (0, 1] 
imply the estimates given in these papers. 



4 Inhomogeneous Problems 

We conclude the discussion of solution representations and a-priori estimates 
with remarks to the inhomogeneous Cauchy problem 

Dv + ^v t = f, v(0,-)=v u v t {0,-) = v 2 (4.1) 
to data v 1 ,v 2 G S(R n ) and a right-hand side f(t,x) G C°°(R + , S(R n )). 



4.1 Solution Representation, Duhamel's Principle 



If we have a nontrivial right-hand side in ( [4,l|) we need entries of the full 
fundamental matrix <!>(£, to,£) to represent the solution. 

To obtain the representation we again apply a partial Fourier transform 
to (4.1). This yields 

vtt + j^vt + \t\ 2 v = f (4.2) 
with data v(0, •) = v\ and vt(0, •) = v 2 - If / = we have 

«(t,0 = $i(t,0,0wi (4-3) 

i=l,2 
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from the discussion of Section ^. Now we assume Vi = 0. We claim that in 
this case 

f)(t,0 = f $ 2 (t,r,0/(r,e)dr. (4.4) 
•/ o 

Indeed we have 



d t / $ 2 (t,r,0/(r,0dr 



Jo 

Jo 



by the initial condition <3? 2 (i,i,f;) = and 
d 2 t /$ 2 (f,r,e)/(r,6dr 



(3t$ 2 )(M, £)/(*> + / 9 t 2 $ 2 (t,T,e)/(r,0dr 



Jo 

from (<9 t <l> 2 )(i,i, £) = 1. Hence if we denote P = <9 2 + jr^d* + |£| 2 we have 
from P$ 2 = 



P / * 2 (t,T,0/(T,f)dT 
JO 

= /(U) + / P$ 2 (t,r,0/(r,0dr = /(*,£)■ 
«/ o 

Obviously function (4.4) satisfies zero initial conditions 

= 0, 



a 



$ 2 (*,r,0/(T,0dT 
$ 2 (t,r,0/(r,e)dr 

$ 2 (o,o,o/(t,e) + 



t=0 



t=0 



a^ 2 (t,r,0/(r,0dr 



0. 



t=o 
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Hence combining both representations by linearity we have 

v(t,0= E *i(*,0,£)«i+ / $ 2 (*,r,0/V,e)dr. 
i=l,2 - 70 



(4.5) 



If we set 

Ki(t,t ,x) = (2 7 r)-fj^J$ i (t,t ^)] ) z = 1,2, (4.6) 

Fourier transform in the sense of 5'(IR n ), we have convoplution representa- 
tion 

v(t,x) = Ki(t,0,x) *vi(x) + K 2 (t,0,x) *v 2 {x) + / i4T 2 (t, r, x) * f(r, x)dr. 

Jo 

(4.7) 



4.2 A-priori Estimates 

We sketch only how to obtain a-priori estimates containing a nontrivial right 
hand side /. We restrict ourselves to the simpler case of estimate for the 
Lq-energy. 

Theorem 4.1. Let u € <S(]R n ). Let further p € (1,2] and q such that 
pq = p + q. Then the estimate 



\\(d t ,V)K 2 (t,T)*u\\ q < 



1 + t 

1 + T 



/ioWs for all < t < t with r = n \ ) . 

J — — \p q J 



The key point is to include also the variable r into the phase variables 
and in the definition of the zones. We use 

Z 1 = {(1+t)|£|>K}, 

Z 2 = {(l + t)\S\>K>Q. + T)\Z\}, 

Z 3 = {K>(l + t)\£\}, 



and follow in each part the reasoning of Section 3.1 





z 3 


z 2 ^ 


Zi 






K 
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As consequence we obtain the following estimate of the L g -energy of the 
solution to ([O]) . 

Corollary 4.2. It holds 

\\(d t ,vwt,-)\\ q 



<(i+*r 



pi i 



Pir+1 + \\v2\\L p , r + [\l + rr 5 \\f(T,-)\\ Lptr dr 
Jo 



(4.8) 



for p G (1,2], g with pq = p + q, r = 

1 fl 1 



and 



max 



n 



P 



1 1 



-, —n I 

P Q 



1 1 



1 



Remark 4.1. The estimate is much weaker than the corresponding estimate 
for the free wave equation, which contains (1+t— r) instead of (l+i)/(l+r). 
The reason for this is the more complex structure of our model problem and 
the absence of the invariance under time translations. 



A Appendix — Relation to Known Cases 
A.l The case fj, = 

The multiplier representation for the solution of ( |1.1| ) obtained in Section ^ 
can be simplified in the case fi = 0. In this case we have p = 1/2 and hence 
the Bessel functions occuring in ( |2.21| ) reduce to trigonometric functions by 

Ji(r) = J—smT (A.l) 

2 V 7TT 

c7_i(r) = \ — cosr. (A.2) 
2 V vrr 

The structure of the multiplier simplifies by the addition theorem of trigono- 
metric functions. This yields 

$i(t,io,£) (A.3) 

TT 



-KlV(l + t)(l+t ) 

Vi ((1 + io)|£|)j7i ((1 + t)|ei) + J-i ((1 + *o)iei) i ((1 + Oiei) 

sin ((1 + toM) sin ((1 + i)|£|) + cos ((1 + t )|£|) cos ((1 + 
cos((t-t )|e|)- (A.4) 
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For the second multiplier one obtains in a similar way 



MMo.O = Sin( V )l;l) . (A.5) 

These multipliers representing the solution of the wave equation are well- 
known. 

The difference t — to is natural in this represention because of the invari- 
ance of the differential operator under time translations. This can only be 
valid for fj, = 0. 

For all even values of [i it is possible to represent the occuring Bessel 
functions by trigonometric functions, but only for fi = 2 the representation 
simplifies significantly, for the reason cf. [A.2j . 

In the general case an addition theorem for Bessel functions like that for 
trigonometric ones is not available. 

A. 2 Transformation to a Klein-Gordon Equation with Vari- 
able Mass 

It is possible to transform the equation 

Uv + ^- t v t = (A.6) 
to a Klein-Gordon type equation. We follow [Mat 76] and [RciOO]. We have 



expi-l f:-^-dT |> = (l + i)-f . (A.7) 



2 7n 1+t 



If we set 



we obtain 



(l+t)~2 W (A. 



= ((1 + f)"2 W ) (i - (1 + t)-$Aw + iff ((1 + t)-$w) t 
= ( - f (1 + tyl- 1 ™ + (1 + t)-2 Wt ) t - (1 + t)-S Aw 

+ - + t)-* _1 w + (i + 

= -f K - 1) (1 + t)-%~ 2 w - + t)-*-V + (1 + t)-$wu 

- (1 + f)"2Affi - £ (1 + t)~2- 2 u; + ^(1 + t)~2- l Wt 



(i + 1 n w + + ^-2^ l 



2(3 



After multiplication with (1+t) 2 we obtain for w the Klein-Gordon equation 

aw + ^zEl(l + t)- 2 w = (A.9) 

with mass tending to zero as t tends to infinity. The sign of the mass term 
depends on the size of /i. For \i G (0, 2) the mass is positive and for fj, > 2 
negative. For the special case /i = 2we obtain a wave equation. 

Remark A.l. This gives an interpretation for the change in the behaviour 
of the energy estimate for ^ = 2. In the forthcoming papers we will develop 
this in more detail for more general classes of weak dissipation. 

Hence if fi = 2 known a-priori estimates for the wave equation (cf . | |5tr7C ] 
W7j| ) lead immediately to 



or 



n-l 



I K*, 01 loo < C e (l + t)-— - 1 (\\v 1 \\ w n + . + \\v 2 \\ wr i+,), (A.lOa) 
IKMII2 <C(||t;i||2 + ||«2|| H -i), (A.lOb) 
and for the derivatives 



(«t(t,0,Vt;(i,0)||oo < Ce(l + t)"^~ 1 (||vi|| VKr+ l +e + \\v 2 \\ W n+e), 

(A.lla) 

\(v t {t r ),Vv(t r ))\\2<C0. + t)- 1 (\\v 1 \\ H i + \\v 2 \\2). (A.llb) 



These estimates are related to the estimates from Section p.2| and [Sj. The 
occurence of the parameter e in these Li-L^ estimates is necessary. The 
space W±(M. n ) itself is not embedded in Loo(lR n ). 
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